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Properties of 1,4-trans poly (isoprene) at ambient conditions are determined by simulations on two 
length scales based on two different models: a full-atomistic and a mesoscopic one. The models are 
linked via a mapping scheme such that one mesoscopic bead represents one chemical repeat unit. 
Melts as well as solutions of several chain lengths were investigated and mapped individually to 
the meso-scale. The resulting models are compared to each other. The meso-scale models could be 
simulated over a large variety of chain lengths and time-scales relevant for experimental comparison. 
Concerning static properties, we determined the persistence length of our systems and the scaling 
behavior of the radius of gyration. The latter was compared to experiments and the agreement is 
satisfactory. Furthermore, we find deviations from Rouse dynamics at all chain lengths at ambient 
conditions. 



I. INTRODUCTION 



Poly (isoprene), better known as rubber, is one of the most commonly used polymers. Its natural form, kautschuk, 
is highly abundant in nature and for technological purposes, synthetic variations are easily industrially polymerised 
from isoprene As an elastomer, its rich field of technical applications comprises tires, textiles, or cable coatings [|). 
Still, some of its static as well as dynamic properties are still not understood on a microscopic basis. 

The time- and length-scales important for the investigation of polymeric materials generally cover a wide range - from 
Angstroms and picoseconds to millimeters and eventually hours. There is no hope to capture all important quantities 
of a system with one single model. Thus, computational methods have to be developed to use the information of more 
detailed models as an input for coarser models allowing to investigate larger systems. The most prominent example 
is the step from quantum chemistry (based on Schroedinger's equation) to semi-empirical atomistic models (based on 
Newton's equation of motion). For many modern applications, however, this step is not sufficient. So, the mapping 
of atomistic systems to meso-scale systems (in which monomers, thought of as "super-atoms", are the smallest units 
of the simulation) became the focus of several studies. 

Coarse-graining (CG) approaches have reached considerable attention over the last decades ||. The basic idea 
of all CG schemes is to separate the system into fast and slow variables where the fast degrees of freedom are 
not of primary interest for the question under study. Thus, a model can be compiled containing the slow degrees 
of freedom only. Recently methods have become available to perform coarse graining in a polymer-specific and 
systematic way [0, M, |L [z| . We refer to a systematic procedure if the identity of the polymer is not lost during the 
mapping procedure. Moreover, the purpose of such techniques is to develop meso-scale models for the simulation of 
large systems of specific polymers. To ease the mapping procedure, two automatic optimization schemes have been 
developed by our coworkers and us: firstly a simplex algorithm technique (^? ], and secondly the iterative Boltzmann 
inversion method which directly targets the problem of structure differences with a physically motivated approach [fl0| . 
Similar ideas have been applied to reproduce experimental data in atomistic simulations 

The goal we pursue here is that two models on different length scales should produce the same distribution functions 
on the larger of the two length scales. Then, we can regard them representing the same polymer. Due to its rich 
applications, we take poly (isoprene) as an example. Additionally we point out some differences to recent findings in 
studies of semi-generic bead-spring models incorporating stiffness as the only local characteristics [O, O . 
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II. INVESTIGATIONS WITH ATOMISTIC DETAIL 

As a basis for our mapping procedure wc need atomistic simulations with good accuracy. We performed simulations 
with appropriately designed force-fields reproducing the atomistic structure and the thermodynamics very reliably. 

A. Melt Investigations 

The atomistic melt simulations have been described in detail in reference [ fTi] . Here, we briefly summarize the 
main characteristics. The simulation box contains 100 oligomers of average length 10 monomers. All chains represent 
trans-poly (isoprene) (cf. Figure [l]). We use a self-developed all-atom force- field resulting in 132 interaction sites for 
a 10-mer [Q; every atom is represented by one interaction site (cf. Tables |]-[n|). The simulations lasted for 1 to 
2 ns at ambient conditions (T — 300 K, p = 101.3 kPa), cf. Table |y|. The model is capable of describing local 
reorientation functions in comparison to NMR measurements and reproduces reasonably the melt structure factor of 
poly (isoprene) [[14, 15 1. 



Atoms connected by any bonding potential did not interact by the Lennard- Jones potential. Additionally, the 
following non-bonded interactions were excluded: all within one monomer, and all C— C, C— H and H— H interactions 
between the second half of the carbons of one monomer (atoms C3, C4, and C5) and the first half of its following 
neighbor (Ci and C2) including the hydrogens connected to them. 

Constant temperature and pressure are ensured using Berendsen's method [ [To] with time constants 0.2 ps for 
temperature and 8 ps for pressure, respectively. The pressure coupling using a compressibility of 2 x 10 -7 kPa -1 was 
employed for the three Cartesian directions independently. All simulations were performed using the YASP molecular 
simulation package |l^] with a time-step of 1 fs and a cutoff for the non-bonded interactions at 0.9 nm. Configurations 
were saved every picosecond. No charges were used through the simulations. 

B. Atomistic Simulations of Poly (isoprene) in Solution of Cyclohexane 

For the atomistic simulations of trans-Pi in cyclohexane we used the same force-field as in the melt simulations except 
that only interactions up to 1 — 4 interactions (inclusive) are excluded. The cyclohexane force-field was developed 
using the automatic simplex method [[l~L| and tested in the dynamics of mixtures of cyclohexene and cyclohexane |p~8| . 
For this study we augmented the cyclohexane force-field with torsion potentials along the ring in order to avoid 
indavertent flipping between the chair and boat configuration fL9| . The densities of the solutions increase slightly 
with concentration of polymer (Tab. ^). They are very close to the density of pure cyclohexane without torsions 
(767.7 kg/m 3 ) pi [lsf , which in turn was optimized against experimental data. We did not reoptimize the Lennard- 
Jones parameters. Therefore the pure cyclohexane system has a slightly lower density as the molecules are less flexible 
and less spherical due to the torsions. In order to decide on the quality of the solvent the radius of gyration 



R G = 4 Yimrf)/ V rrn, n = Ik - (At) (1) 



and the end-to-end length are important characteristics. The mean square end-to-end distance v/ (Re-e) ^ s measured 
between the terminal carbons of the chain. Additionally the hydrodynamic radius Rh was calculated as the first 



inverse moment of the distance vectors 20 




Here, N represents the number of monomers of the chain and the distance between two arbitrarily chosen super- 
atoms i and j. 

We do not find a systematic variation of the radius of gyration with concentration. In the system with two chains 
in a solvent of 500 molecules the two chains behave differently whereas in the higher concentrated systems the radii 
of gyration are very close although they were also started with the same initial configurations as the other two chain 
system. The auto-correlation functions of the gyration radii, however, decays on the scale of less than a hundred 
pico-seconds. Thus, the systems are equilibrated even though the error bars are quite large. In order to decrease them 
substantially much larger simulations would be necessary. This was not the purpose of the atomistic simulations as 
the long time-scales are left to the coarse-grained model to be developed here. 
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In Figure g we show the decay of the reorientation correlation function of the end-to-end vector of the whole chains 
in solution. We see that this function overall follows an exponential tendency and the correlation time is in the order 
of 2 — 4 ns. This supports once again the assumption of equilibration. Additionally the overall reorientation becomes 
slower with increasing concentration. The dynamics of the solvent is found to slow down and become more anisotropic 
with increasing concentration |1S|]. 



III. THE MAPPING PROCEDURE 



The mapping procedure is the same for both the melt and the solution system. It has already been described in full 
detail in reference []l0|| , so we limit ourselves to a short summary. Distribution functions from atomistic simulations 
are the basis of the mapping, since the mesoscale model is optimized against them. In this study, Ml (for the melt) 
and the combination of SI and S3 (for the solution) were chosen as parent atomistic simulations. The mapping itself 
is illustrated in Figure |l|. Each chemical repeat unit is replaced by one super-atom (interaction center), located in 
the middle of the atomistic bond between successive chemical repeat units. Consequently an atomistic A-mer will be 
coarse grained to a A — 1-mer. This is meaningful because the such mapped super-atoms are much more spherical 
compared to the case in which one would have chosen the center of mass of the real monomer. The advantage is that 
the mapped chains generate well distinguishable peaks for the various intramolecular distributions. Especially when 
applying the model to longer chains, the error introduced due to the missing rests are negligible. 

In this contribution, we use exclusively the pressure-corrected CG force field which has been optimized by the 
Inverted Boltzmann method. It is fully described in Reference JToj ] . The idea, however, shall shortly be illustrated 
here. The Inverted Boltzmann method utilizes the differences in the potentials of mean force between the distribution 
functions generated from a guessed potential and the true distribution functions to improve the effective potential 
successively. Imagine we would like to derive an effective non-bonded potential Voo(r) from a given tabulated start 
potential Vo(r), targeting to match the radial distribution function g^r). Simulating our system with Vo(r) will 
yield a corresponding go(r) which is different from <?oo(t"). The potential needs to be improved, which is done by a 

correction term — fc^Tln ( 3 g ff ) ■ This step can be iterated as follows by: 

V j+1 (r) = Vj{r)-k B Tin(^j (3) 

until 

/target = / uj(r) (goc(r) - 9j{r)f dr . (4) 

falls below an initially specified threshold. (We apply w(r) = exp(-r) as weighting function in order to additionally 
penalize deviations at small distances.) Pressure corrections can be introduced in polymer systems by adding a week 
linear potential term to the attractive outer part of Vj(r). One can then post-optimize the structure of the system 
according to the above scheme in turn with occasional linear additions until the pressure matches the one from the 
atomistic system, too. 



IV. MESOSCALE SIMULATIONS 



A. Technical Simulation Details 



Both Brownian Dynamics (BD, for the melt) and Monte Carlo (MC, for the solution) programs have been applied. 
All BD runs were performed in the NVT ensemble, with values for density and temperature corresponding to the 
values of the parent atomistic simulation. The systems consisted of an orthorhombic box employing periodic boundary 
conditions. Langevin equations of motion were integrated by the velocity Verlet algorithm with a time step At = 
O.Olr [2l| | at a dimensionless temperature of T* = 1. This temperature was maintained by the Langevin thermostat 
with friction constant T = 0.5t _1 |^2). For the MC simulations representing the solution case, a single chain program 
was used p3[ . Here, 10 5 accepted warm-up moves were carried out before a production run of 10 6 accepted moves 
was started. For analysis, every 500th configuration was stored, both in BD and MC simulations. 
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B. Static properties in the melt and in solution 



The results of the CG simulations are summarized in Table [Vj| in case of the poly (isoprene) melts and in Table VII 



in case of the solutions. For both, the radius of gyration Rq and the hydrodynamic radius Rh is correctly reproduced 
compared to the values of the parent atomistic simulations (within standard errors). 

Extensive Brownian Dynamics simulations were carried out to investigate the static and the dynamic behavior of 
poly (isoprene) melts. The coincidence (within statistical fluctuations) of the values for D center and D com , representing 
the diffusion coefficient of the central monomer of the chain and their center of mass, respectively, shows that the 



system is fully equilibrated. Let us here concentrate on the static properties. In Table VI we list the values for the 
radius of gyration Rq, the hydrodynamic radius Rh, their ratio and the ratio R\jR G which provide information 
about the extension of the chains. The scaling behavior of Rq is shown in Figure [| A fit for N > 35 yields a scaling 
exponent of v = 0.53, which is slightly larger than the ideal value for melts of ve = 0.5, suggesting that the chains are 
not sufficiently close the the limit of infinitely long chains. This is also supported by the ratio Rq/Rh which is well 
below the infinite limit of around 1.5. The fact that R\jR G — 6.0 holds for all chain lengths, however, shows that 
the static behavior corresponds to random walks. 



Table VII holds the results of the dissolved poly (isoprene) systems. Since we applied MC, static properties are 
listed exclusively. The scaling behavior of Rq is also presented in Figure [|. With a fit for the same region as in the 
melt case (N > 35) we obtain a scaling exponent of v = 0.60. This is only slightly larger than the theoretical value 
v K s = 0.588 for polymers in good solvents due to the finite chain length. The ratio R^/Rq decreases monotonically 
with chain length. This indicates that the intramolecular interactions of poly (isoprene) in the trans state locally 
stretch the chains. 

To the best of our knowledge only chain size measurements of cis-poly (isoprene) in solution of cyclohexane have 
been performed |24|, |25| . Our data coincides within the order of magnitude with the extrapolated experimental data 
of Tsunashima et al. to our longest chains. Applying their relation for the gyration radius we calculate 2.76 nm for 
a chain length of 100 monomers in O-solution. This is almost exactly our value for the chains of this length in the 
melt where also random walk statistics is supposed to apply. Davidson et al. measured the radius of gyration of a 
mixture of cis and trans poly (isoprene) under good solvent conditions by various techniques including wide angle 
light scattering for molecular weight range of 112,000 g/mol and higher. p4j This corresponds to chain lengths of 1650 
monomers and longer. Extrapolating the simulation data to the two shortest chains experimentally available we get 
23.7 nm for 112,000 g/mol compared to 17.7 nm experimentally and 26.5 nm for 164,000 g/mol against 19.4 nm 
experimentally. Thus, we overestimate the gyration radii by about 35%. However, one has to keep in mind that the 
experimentally investigated system contains a mixture of the trans and cis conformer. The experimentally observed 
scaling exponent is 0.545 and therefore lower than our exponent of 0.6 and also lower than the theoretically predicted 
exponent of 0.588. 

Next, the persistence length ^ p of poly (isoprene) was estimated by the decay of the orientation correlation along 
the chain backbone. The point where this correlation arrives at y = exp(— 1) is a good approximation of Z p 

(uk(n + i) u k (n)\ = exp(-^M . (5) 



\ v J ' v 7 VI. 

Here, lb is the bond length and Ukin) represents the unit vector along the chain, centered at monomer n and pointing 
to monomer k. The result is shown in Figure ^. The bonds are completely decorrelated beyond the forth neighbor for 
the melt. For the persistence length, l p sa 1.25Zb(= 0.59 nm) can be estimated. As expected, this is independent of 
chain lengths, as all curves fall on top of each other. The same holds for the solution case. Here the persistence length 
could be estimated to be l p w 2.7Zb(=1.27 nm), i.e. more than twice as large as for the melt. So, the concentration 
of poly (isoprene) chains has a strong influence on the intrachain statistics. The dense packing in the melt favors a 
quick decorrelation between neighboring chain monomers. 

Finally, we can compare some of the above results to data for a melt of generic bead-spring chains with bond-angle 
stiffness |la. It has already been used to successfully map the dynamic behavior of atomistic poly (isoprene) the 
meso-scale ||13|], even though it just comprised one simple bond angle potential in addition to a FENE bead spring 
interaction for the intramolecular part: 

Vanglc = X ■ [1 + Uj ■ Uj+l] k B T . (6) 

The vector u corresponds to a normalized bond vector and x is the force constant to determine the stiffness. The 
scalar product corresponds to the cos ine o f the angle along the backbone. The comparison of static properties of 



the two models is presented in Table VIII . One observes that most properties yield similar values, in case of Rq, 
Rl/R G , and ^ p even within statistical error. However, the end-to-end distances deviate more strongly, showing that 
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the realistic poly (isoprene) chains are more coiled than generic chains. This is also reflected in the Kuhn segment 
length Zkj which defines the distance beyond which one expects universal chain behavior: 

i (jje ) (jje ) f-i j irj\ 

~ R[ max) ~ {N-l)h ~ °° b ' [ ' 

It is significantly shorter for the realistic model compared to the generic one. A reason for that could be the missing 
torsional potential in the latter. Also, the dynamical behavior turns out to differ significantly, as will be shown in the 
next section. 

Although the mesoscopic model is optimized against the atomistic structure, some structural differences remain. 
Comparing the pair distribution functions of centers-of-mass of whole chains (cf. Figure |^) we see that the mesoscopic 
model chains show less structure on short distances. The mesoscopic 9-mers show just a correlation hole whereas the 
atomistic 10-mcrs exhibit a structured correlation hole. Moreover we see the strong difference between a mesoscopic 
9-mer and a 50-mer. In the latter case the correlation hole is much less pronounced. This can be understood as the 
bigger a chain becomes the less volume of its ellipsoid of gyration it actually occupies and the more overlap of such 
ellipsoids is possible. 



C. Dynamical Properties in the melt 

The Rouse model |^6| is widely accepted for the dynamics of short flexible chains in the melt. The assumptions it is 
based on are the screening of excluded volume and that all local intra-chain interactions can be mapped onto a Kuhn 
segment length [f27| . Poly (isoprene) has been experimentally tested against the Rouse model Q and the agreement 
was not too good at room temperature. Simulations at the atomistic model found for the oligomers of length 10 a 
reasonable description by the Rouse model at the elevated temperature of T = 413 K fl3[ . 

For long chains deviations from the Rouse model are expected due to mutual topological constraining of chains. 
Such deviations could then be attributed to entanglements since recently it has been shown that even a weak presence 
of stiffness can dramatically decrease the onset of entanglement influences |l2| . 

Two observables are commonly used to decide on the validity of the Rouse model: the mean-squared displacements, 
and the so-called Rouse modes. Figure || shows the g\ function which is defined as 

gi(t) = ((f l (t + t )~r l (t )) 2 ) to , (8) 

where fi(t) denotes the position of the central monomer of a chain at a given time. The central monomer is chosen to 
suppress end-effects as good as possible. For short times (t < tr) one expects all curves to fall on top of each other; 
tr is the Rouse time with all internal degrees of freedom are relaxed (below). At such short times the individual 
monomer does not yet have to drag the whole chain with it, but an increasing part. At t rj Tr(N) the different 
curves leave the master curve as the increasing neighborhood has reached the chain length. For longer chains Figure ^ 
shows clearly such a Rouse regime where gi(t) cx t - 5 . Additionally, for all chain lengths N < 85 we find the free 
diffusion limit gi(t) cx t reached and therefore the systems equilibrated. This shows that the local dynamics is at least 
Rouse- like with the connectivity leading to a time dependent diffusion coefficient for t < tr. This in turn leads to 
the sub-diffusive behavior observed. For longer chains the local regime lasts longer which clearly hints towards the 
connectivity behind the sub-diffusivity. 

The actual Rouse model relies on the independence of the Fourier modes of the chains. We compare the decay of 
these Rouse modes X p 

i=0 v ' 

According to theory the curves for the different modes at a given chain length should collapse if the time and the 
correlation function are rescaled by the squared mode index p0| . For intermediate chain lengths we find this scaling 
to work reasonably well. For longer chain lengths, however, the higher modes are apparently faster than predicted 
by the Rouse model. This is first hint towards entanglements. It has been observed in simulations of other polymer 
models that especially the lower modes are likely to be the first to be influenced by entanglements O, M . 



Table IX shows the Rouse time tr derived by an exponential fit to the first Rouse-Mode X\. We clearly sec finite 
chain length effects. For longer chains N > 100 these effects start to become weaker. Figure ^ compares the scaling 
for several chain lengths. This suggests that poly (isoprene) melts at room temperature do not completely behave 
according to Rouse theory. 
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We are not yet able to reach into the entanglement regime. Especially for the longest chains the scaling of the Rouse 
modes for different chain lengths onto each other apparently holds. It is interesting to note that we see finite chain 
length effects for chains of moderate length. Figure [7]d is a way to tell oligomers from polymers as the self-similarity 
is obeyed in the polymer case. For oligomers the Rouse friction is still chain length dependent. This was not found 
in a recent study of model chains only incorporating stiffness. In that case the Rouse model was obeyed almost from 
the very beginning (N > 15). This is another important finding in the view of modeling long oligomers or short 
polymers. For chains with chemical identity we find a regime where the chain nature is already important but the 
generic polymer models are not yet valid. 

V. CONCLUSIONS 

Atomistic simulations of trans-poly (isoprenc) in the melt and in solution have been simulated and independently 
coarse-grained to models with one interaction center per monomer. The local scale models arc verified against 
experimental data where available. In the melt a slowdown of the dynamics is found as expected. For the solution 
we were able to equilibrate oligomers of length 15 monomers until the end-to-end reorientation correlation function 
decayed to less than about 30% of its initial value. Note, that we did not need heavy hardware resources to generate 
the statistical precision we present here. 

The two independently coarse-grained models for solution and melt respectively are distinctly different. It is 
impossible to represent a chain in such strongly different environments on the coarse grained level with a generic 
potential. This is especially true as we aimed to get rid of the solvent completely in the solution case and melt and 
solution are very well known to behave strongly different in the long chain limit. 

The coarse-grained models could be compared against well-known long polymer scaling theory. Additionally, com- 
parison to experiments yielded a difference of 35%. The Rouse behavior is for this model recovered for chains of more 
than around 80 monomers. In summary, we find that the proposed coarse graining procedure allows us to produce 
models which are reliably reproducing properties of real polymers. We were not yet able to reach into the reptation 
regime. 
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FIG. 1: Illustration of the mapping of trans- 1,4-poly (isoprene) from the atomistic to the mesoscopic level. Each chemical 
repeat unit is represented by one super-atom. As center of these super-atoms, we choose the middle of the atomistic bond 
between two successive chemical rpnpnt nnit« Tl-n« i« lla pfiil hppnnsp tViP rpsnltino- mnrmpH rtiains o-pnerate well distinguishable 
peaks for the various intramolec 




time t [ns] 



FIG. 2: Reorientation correlation function of whole chains in solution in semilogarithmic representation. We show the second 
Legendre polynomial P2 of the reorientation angle. The shortcuts SI, S3, and S4 are defined in Table 



9 



[nm] 


_ . _ ___ 60 

RprirpQQinn* v— D ?7Hx 

i ity i uojiui i . y — \j . l— i ua 




/ration 


Regression: y=0.245x 053 22 




CD 






Radius of 


• 


#CG melt 
OCG solution 




10 


100 



Chain length [monomers] 

FIG. 3: Scaling behavior of CG poly (isoprene) chains in melt (filled circles) and solution (open diamonds). 



10 



"D 

CD 
N 



O 

c 



CD 

i— 
o 
o 
c 
o 



c 

CD 



0.8 



0.6 



« 0.4 - 



0.2 



-0.2 



* ' 1 

w 

■ \\ X = 1 

-\\ 
\^ 

\ 

\ 


.25 

\ 

■\ 
\ 
■-\ 


x = 2.7 

G> o CG 35-mer 

A- -A CG 85-mer 

y = exp(-1) 




\ 

\ 

\ 

\ 

\ 

\ 

k. 


solution "-^-^j 

■ — 6 —4i 






melt 



Distance [# monomers] 



FIG. 4: Approximation of the persistence length l p of some arbitrarily picked CG poly (isoprene) systems. Using eq. |f| the 
decay of the orientation correlation function to y = exp( — 1) approximates l p . In the melt, this yields / p « 1.25Zb and for the 
2.7l b . 



solution, l p 



1.25 




Radius r[nm] 



FIG. 5: Pair distribution function g(r) for different chain models. We show two different atomistic 10-mers (Ml and M2), a 
mesoscopic 9-mer, and a mesoscopic 50-mer. 
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FIG. 6: The gi-function, mean-squared displacement of central monomers for the melt simulation at chain lengths N = 
9, 15, 20, 35, 40, 45, 50, 85, 120 from top to bottom. As guides to the eye we show lines indicating i 0,5 and t 1 . 
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FIG. 7: Decay of the autocorrelation functions of the Rouse Modes for different chainlengths in Rouse scaling. Figure a) N = 35 
b) N = 65, c) N = 85. The symbols are described in figure b. d) Comparison of the decay of the autocorrelation function for 
different chain lengths. 



TABLE I: Angles (equilibrium values <j>o and potential stren gth k) and bond lengths li for atomistic melt simulations of poly 
(isoprene). The bond lengths are constrained using SHAKE |30|. 





angle 




</>o 


fc[kJ/(mol*rad 2 )] 


bond 


h H 


Ci 


-Ca- 


c 3 


128.7 


250 




Ci — C 2 


0.150 


c 2 


-C 3 - 


c 4 


124.4 


250 




C 2 = C 3 


0.1338 


C 2 


-C 3 - 


C 5 


120.2 


250 




C3 — C4 


0.151 


c 4 


-C s - 


C 5 


115.4 


250 




C3 — Cs 


0.1515 


c 3 


-Cs - 


Ci 


114.5 


250 




Cs — Ci 


0.155 


c 5 


-Cr- 


c 2 


112.7 


250 




C-H 


0.109 


c- 


C ap 3 


-H 


109.5 


250 








Ci 


-C a - 


-H 


114.4 


250 








c 3 


-C a - 


-H 


114.4 


250 








H 


-C- 


H 


109.5 


250 
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TABLE II: Force-field parameters of the non-bonded interactions for atomistic melt simulations of poly (isoprene). m is the 
atom mass, a the interaction radius, and e the interaction strength. 



atom 


m [amu] 


a [nm] 


epsJ/mol] 


C sp2 


12.01 


0.321 


0.313 


C sp3 


12.01 


0.311 


0.313 


H 


1.00782 


0.24 


0.2189 
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TABLE III: Force-field parameters of the dihedral angles for atomistic melt simulations of poly (isoprene). The identifying 
numbers are found in figure Fl 



torsion equilibrium an; 


de barrier 


periodicity 


ro [degrees] 


k T [kj / mol] 


i 


1 


5.2 


1 


1 


-7.4 


2 


1 


10.0 


3 


2 180 


9.7 


1 


2 180 


14.1 


3 


3 


-21.1 


1 


3 


-12.3 


2 


3 


0.5 


3 



TABLE IV: Summary of the characteristics of the atomistic melt simulations. 



system T[K] 


isim[ps] 


M„ 

M„ 




Ml 300 


1184 


1.00 


890 


M2 300 


2012 


1.05 


917.4 


M3 300 


1737 


1.05 


916.8 


M3 413 


792 


1.05 


826 



TABLE V: Thermodynamic and static properties of the solution. Details of the different systems and the numbers with which 
they are referenced in the following. Np is the number of oligomers (C75H122) and Nc the number of cyclohexane molecules, 
c is the concentration in weight % polymer. i s i m is the simulated time for the systems. 



system Np 


N c c 


P [kg/m J ] 


U\ m [ns] 


Rh h 


Re [nm] 


Rti 1 Rg 


SO 


500 0.0% 


756.3 


1.0 








SI 1 


250 4.6% 


764.2 


11.25 


1.33 ± 0.12 


1.21 ± 0.20 


6.1 


S2 1 


500 2.4% 


757.5 


5.58 


1.40 ± 0.05 


1.26 ± 0.10 


6.0 


S3 2 


500 4.6% 


762.5 


7.81 


1.34 ± 0.10 


1.23 ± 0.20 


6.5 


S4 2 


250 8.9% 


768.2 


11.56 


1.37 ± 0.09 


1.29 ± 0.15 


7.2 



TABLE VI: Data of the mesoscopic Brownian Dynamics poly (isoprene) melt simulations. The pressure is stated in reduced 
units | pl[ . All systems are simulated at a temperature of T = 300 K and at a density of p v — 7.25 monomers/nm 3 and 
correspond to the atomistic simulation Ml. The static properties come with a standard error of w 5%, the dynamic properties 
with an error of « 15%. 



chain length N 


9 


15 


20 


35 


40 


45 


50 


65 


85 


100 


120 


number of chains 7V C 


100 


150 


150 


100 


150 


150 


80 


130 


50 


120 


100 


pressure p* 


±0.0 


-0.01 


-0.04 


-0.15 


-0.07 


-0.08 


-0.08 


-0.08 


-0.09 


-0.09 


-0.09 


simulation time 


3.0 


4.0 


4.0 


7.2 


6.0 


6.0 


8.0 


6.0 


8.0 


6.0 


6.0 


[10 6 time steps] 
























Rg [nm] 


0.75 


1.03 


1.20 


1.63 


1.75 


1.87 


1.97 


2.25 


2.63 


2.85 


3.14 


Rh [nm] 


1.11 


1.19 


1.27 


1.50 


1.60 


1.65 


1.69 


1.98 


1.98 


2.21 


2.42 


Rg/Rh 
Re/ Rc 


0.68 


0.86 


0.94 


1.09 


1.10 


1.13 


1.17 


1.14 


1.33 


1.29 


1.30 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


-Dcontrc [10~ 6 Cm 2 /s] 


16.6 


10.7 


6.7 


4.2 


2.5 


2.3 


2.9 


1.3 


1.3 


0.8 


0.6 


D COI n [10- 6 cm 2 /s] 


16.0 


10.9 


6.7 


4.1 


2.5 


2.3 


2.8 


1.3 


1.2 


0.8 


0.6 
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TABLE VII: Static properties (with standard errors of ~ 5%) as obtained from Monte-Carlo simulations of mesoscopic poly 
(isoprene) solutions. The system density and temperature corresponds to the parent atomistic systems SI and S3. 



chain length N 


9 


15 


35 


50 


85 


100 


120 


Rg [nm] 


0.91 


1.33 


2.37 


2.92 


4.11 


4.48 


4.98 


Rh [nm] 


1.23 


1.40 


1.98 


2.32 


3.06 


3.31 


3.63 


Rg/Rh 

Re/Rh 


0.74 


0.95 


1.20 


1.26 


1.34 


1.35 


1.37 


7.3 


7.1 


6.7 


6.5 


6.5 


6.4 


6.4 



TABLE VIII: Results of some static properties for a realistic meso-scale model of poly (isoprene) in a melt and a more 
generic bead-spring model with bond-angle stiffness. Rg is the radius of gyration, R e the end-to-end distance of the chain, l p 
corresponds to the persistence length and Ik to the Kuhn segment length. 



Melt Model 


[nm] 


[nm 2 ] 


rI/Rq 


[monomers] 


Ik 
[nm] 


x = 1.5, N = 20 P] 
CG PI, N = 20 


1.28 
1.20 


10.3 
8.7 


6.3 
6.0 


1.2 
1.25 


1.2 
1.0 


x = 1.5, N = 50 Q 
CG PI, N = 50 


2.20 
1.97 


29.5 
23.3 


6.2 
6.0 


1.3 
1.25 


1.3 
1.0 



TABLE IX: Fits of the Rouse time using the decay of the first Rouse mode for different chain lengths. 



N 


TR 


tr/N 2 


35 


743 


0.60 


50 


1946 


0.78 


65 


3722 


0.88 


85 


9300 


1.29 


100 


13800 


1.38 


120 


17700 


1.23 



1.25 




Radius r [nm] 



2.5 



2 - 



1.5 - 



LL 
Q 

1 



0.5 - 



--- Full RDF 

RDF-a (1 neighbour excluded) 

RDF-b (2 neighbours excluded) 




0.25 



• ar m «* 



0.5 



0.75 1 1.25 
Radius r [nm] 



1.5 



T3 
CD 
N 

"ro 
E 

O 

c 

(£ 

3- 
°l 

a 



0.02 



=- 0.015 



0.01 
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CD 
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